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Abstract 

An upper bound on the error probability of specific lattices, based on their distance-spectrum, is constructed. 
The derivation is accomplished using a simple alternative to the Minkowski-Hlawka mean-value theorem of the 
geometry of numbers. In many ways, the new bound greatly resembles the Shulman-Feder bound for linear codes. 
Based on the new bound, an error-exponent is derived for specific lattice sequences (of increasing dimension) over 
the AWGN channel. Measuring the sequence's gap to capacity, using the new exponent, is demonstrated. 



I. Introduction 

(N- 

For continuous channels, Infinite Constellation (IC) codes are the natural coded-modulation scheme. The 
c3 encoding operation is a simple one-to-one mapping from the information messages to the IC codewords. 
The decoding operation is equivalent to finding the "closest'^ IC codeword to the corresponding channel 
output. It has been shown that codes based on linear ICs (a.k.a. Lattices) can achieve optimal error 
performance HI, ED. A widely accepted framework for lattice codes' error analysis is commonly referred 
to as Poltyrev's setting In Poltyrev's setting the code's shaping region, defined as the finite subset of 
\ the otherwise infinite set of lattice points, is ignored, and the lattice structure is analyzed for its coding 
(soft packing) properties only. Consequently, the usual rate variable R is infinite and replaced by the 
Normalized Log Density (NLD), 5. The lattice analogous to Gallager's random-coding error-exponent 
BH, over a random linear codes ensemble, is Poltyrev's error-exponent over a set of lattices of constant 
CSj \ density. 

Both Gallager and Poltyrev's error-exponents are asymptotic upper bounds for the exponential behavior 
of the average error probability over their respective ensembles. Since it is an average property, examining 
a specific code from the ensemble using these exponents is not possible. Various upper error bounding 
O ; techniques and bounds for specific codes and code families have been constructed for linear codes over 
C^l ■ discrete channels 0, while only a few have been devised for specific lattices over particular continuous 
channels J6j. 

The Shulman-Feder Bound (SFB) [7] is a simple, yet useful upper bound for the error probability of a 
^ ■ specific linear code over a discrete-memoryless channel. In its exponential form it states that the average 
error probability of a g-ary code C is upper bounded by 

p e (C) < e -M*+^) CD 

A^c(r) e nR 

a = max . . . , „ (2) 

r K{t) e nR - 1 

where n, R, and E r (R) are the dimension, rate, and random-coding exponent respectively, and Nc{t) 
and M r {r) are the number of codewords of type r for C and for an average random code, respectively 
(i.e. distance- spectrum). The SFB and its extensions have lead to significant results in coding theory. 
Amongst those is the error analysis for Maximum Likelihood (ML) decoding of LDPC codes [H. The 
main motivation of this paper is to find the SFB analogue for lattices. As such it should be an expression 
that upper bounds the error probability of a specific lattice, depend on the lattice's distance spectrum, and 
resemble the lattice random-coding bound. 

The main result of this paper is a simple upper bounding technique for the error probability of a specific 
lattice code or code family, as a function of its distance- spectrum. The bound is constructed by replacing 
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the well-known Minkowski-Hlawka theorem [9] with a non-random alternative. An interesting outcome 
of the main result is an error-exponent for specific lattice sequences. A secondary result of this paper 
is a tight distance- spectrum based, upper bound for the error probability of a specific lattice of finite 
dimension. 

The paper is organized as follows: Sections II and III present the derivation of the Minkowski-Hlawka 
non-random alternatives, section IV outlines a well-known general ML decoding upper bound, section V 
applies the new techniques to the general bound of section IV, and section VI presents a new error-exponent 
for specific lattice sequences over the AWGN channel. 



II. Deterministic Minkowski-Hlawka-Siegel 

Recall that a lattice A is a discrete n-dimensional subgroup of the Euclidean space R n that is an 
Abelian group under addition. A generating matrix G of A is an n x n matrix with real valued coefficients 
constructed by concatenation of a properly chosen set of n linearly independent vectors of A. The 
generating matrix G defines the lattice AbyA = {A:A = C7u, u e Z"). A fundamental parallelepiped 
of A, associated with G is the set of all points p = Y^7=i Ui 9i where < Ui < 1 and {gi}™ =1 are the basis 
vectors of G. The lattice determinant, defined as det A = |degC7|, is also the volume of the fundamental 
parallelepiped. Denote by (3 and 5 the density and NLD of A respectively; thus (3 = e nS = (det A) -1 . 

The lattice-dual of the random linear codes ensemble, in finite-alphabet codes, is a set of lattices 
originally defined by Siegel ifTOl . [fTTTl . for use in proving what he called the Mean- Value Theorem (MVT). 
This theorem, often referred to as the Minkowski-Hlawka-Siegel (MHS) theorem, is a central constituent 
in upper error bounds on lattices. The theorem states that for any dimension n > 2, and any bounded 
Riemann-integrable function g(X) there exists a lattice A of density f3 for which 



AeA\{0} 



g(\) < [ g(^)dx = p f p(x)dx. (3) 

J R n \e J J R n 



Siegel proved the theorem by averaging over a fundamental sell of all n-dimensional lattices of unit 
density. The disadvantages of Siegel's theorem are similar to the disadvantages of the random-coding 
theorem. Since the theorem is an average property of the ensemble, it can be argued that there exists 
at least a single specific lattice, from the ensemble, that obeys it; though finding that lattice cannot be 
aided by the theorem. Neither can the theorem aid in analysis of any specific lattice. Alternatives to ©, 
constructed for specific lattices, based on their distance- spectrum, are introduced later in this section. 

We begin with a few definitions, before stating our central lemma. The lattice A always refers to 
a specific known n-dimensional lattice of density (3, rather than A which refers to some unknown, yet 
existing n-dimensional lattice. The lattice A is the normalized version of A (i.e. det(A ) = 1). Define 
the distance series of A as the ordered series of its unique norms {\j}JL , such that Ai is its minimal 
norm and A — 0. is defined for A respectively. The normalized continuous distance-spectrum 

of A is defined as 

oo 

N(x) =Y^Mj8{x- %) (4) 
j'=i 

where {J\Ij}°? =1 is the ordinary distance- spectrum of A , and 5(-) is the Dirac delta function. Let T denote 
the groupj of all orthogonal n x n matrices with determinant +1 and let /i(7) denote its normalized 
measure so that j T d^(^j) = 1. The notation 7A0 is used to describe the lattice generated by 7G, where 
G is a generating matrix of the lattice A . 

2 Let T denote the multiplicative group of all non-singular n x n matrices with determinant 1 and let <!> denote the subgroup of integral 
matrices in T. Siegel's fundamental set is defined as the set of lattices whose generating matrices form a fundamental domain of T with 
regards to right multiplication by $ (see section 19.3 of (9))- 

3 This group, consisting only of rotation matrices, is usually called the special orthogonal group. 
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Our central lemma essentially expresses Siegel's mean-value theorem for a degenerate ensemble con- 
sisting of a specific known lattice A and all its possible rotations around the origin. 

Lemma 1. Let Ao be a specific n- dimensional lattice with NLD 5, and g{\) be a Riemann-integrable 
function, then there exists an orthogonal rotation 7 such that 



E aw < [ m*\\)g 

a \ mi ->R n 



with 



Ae 7 A \{0} 



N(x) 
nVnX 71 - 1 





: x > 
: x < 



(5) 



(6) 



where V n is the volume of an n-dimensional unit sphere, and ||-|| denotes the Euclidean norm. 

Proof: Let denote the subspace of all points 9 in the n-dimensional space with ||0|| = 1, so that 
6 is the surface of the unit sphere. Let fi(6) denote the ordinary solid-angle measure on this surface, 
normalized so that J" e <i/i(0) = 1. We continue with the following set of equalities 



J £ g{X)dp{ 1 ) 



Ae 7 A \{0} 



AeA o \{0} • 



dn{0) 




dn(9)dR 



dfi(0)dV n R n 



R"\{o} nV n || 
K(\\x\\)9 



N(\\*\\)_ n 

n-l 9 



X ' 

— r ) C?X 

e 5 



R" 



<ix. 



(7) 



where the third equality follows from the definition of T and and the measures 71(7) and fi(0), the 
fourth equality is due to the circular- symmetry of the integrand, and the sixth equality is a transformation 
from generalized spherical polar coordinates to the cartesian system (see Lemma 2 of ifTTTO . 

Finally there exists at least one rotation 7 G T for which the sum over 7A0 is upper bounded by the 
average. ■ 

The corollary presented below is a restricted version of lemma Q] constrained to the case where the 
function g(X) is circularly-symmetric, (i.e. g(X) = g(||A||)). To simplify the presentation, it is implicitly 
assumed that g(A) is circularly-symmetric for the remainder of this paper. It should be noted that all 
results presented hereafter apply also to a non-symmetric g(\) with an appropriately selected rotation 7 
of A . 



Corollary 1. Let A be a specific n-dimensional lattice with NLD 5, and g{\) be a circularly-symmetric 
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Riemann-integrable function, then 

0(A) = / 9T(||x||)s(*W (8) 



ASA \{0} 

with 9T(x) as defined in lemma [7] 

Proof: When (7(A) is circularly-symmetric, 



/ £ ^(A)dAi(7) = E / gfrWtii) 

"' T Ae 7 A \{0} AeA \{0}" /r 

E gWf dtii) 



AeA \{0} 

= E 2(A) (9) 

AGA \{0} 

■ 

The right-hand side of © can be trivially upper bounded by replacing %l(x) with a suitably chosen 
function a(x), so that 

J m(\\x\\)g (J) rfx < jf ^ a(||x||)p dx. (10) 
Provided the substitution, it is possible to define the following upper bounds: 

Theorem 1 (Deterministic Minkowski-Hlawka-Siegel (DMHS)). Let A be a specific n-dimensional lattice 
of density B, g(X) be a bounded Riemann-integrable circularly -symmetric function, and at(x) be defined 
such that (flOl) 15 satisfied, then 



y g(A) < B max a(x) / g(x)<ix (11) 
AeA \{o} b><^ J J** 

where A max is the maximal ||x|| for which g(x) ^ 0. 

Proof: Substitute a specific a(x) for %l(x) in ©, and upper bound by taking the maximum value of 
a(x) over the integrated region, outside the integral. ■ 

Theorem 2 (extended DMHS (eDMHS)). Let A be a specific n-dimensional lattice with density 8, 
g(X) be a bounded Riemann-integrable circularly -symmetric function, a(x) be defined such that (flOl) is 
satisfied, and M be an positive integer then 

> q(X) < 8 min > maxa(x) / o(x)dx 1 (12) 



AeA /{0} 1 JJ ' =1 \j=i 



with 



Hj = {x :x> 0, e*Rj-i <x< e 5 R 3 } (13) 
TZj = {x : x e E n , < ||x|| < i^} (14) 



where {i? J }yl 1 w an ordered set of real numbers with R = an J -Raj = A max , where A max w f/ze maximal 
||x|| /or which p(x) 7^ 0. 

Proof: Substitute a specific a (a;) for 0T(a;) in ([8]) and break up the integral over a non-overlapping set 
of spherical shells whose union equals R". Upper bound each shell integral by taking the maximum value 
of a(x) over it, outside the integral. Finally, the set of shells, or rather shell-defining radii is optimized 
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such that the bound is tightest. ■ 
The eDMHS may be viewed as a generalization of DMHS since for M = 1 it defaults to it. In addition, 
when M — > oo the eDMHS tends to the original integral. Clearly, both bounds are sensitive to choice of 
a{x), though one should note, that for the same choice of a(x) the second bound is always tighter. 

The bounds shown above are general up to choice of a(x), and clarify the motivation for the substitution 
in (fTQb . Careful construction of a(x) along with the selected bound can provide a tradeoff between tightness 
and complexity. The next section presents a few simple methods for construction of the function a(x), 
and their consequences. 

III. Construction of a(x) 

Maximization of the right-hand-side of ©, by taking out the maximum value of outside the 

integral, is not well-defined. This is since 9T(x) is an impulse train. The motivation of this section is to 
find a replacement for 9T(x) that enables using this maximization technique whilst retaining a meaningful 
bound. 

Let's assume that g(X) is monotonically non-increasingF| in ||A||, and define N'(x) to be a smoothed 
version of the normalized continuous distance-spectrum, selected such that it satisfies 

rr rr 

/ N(R)dR< N'(R)dR Vr G (0, e 5 A max ]. (15) 
Given the above, a(x) can be defined by expanding (flOl) as follows: 

e s \ n 



1 ot(h)9 ■>*>»(?)«* 

r 4Amax fR\ 
-L N ^{7) dR 



[ «(IMI)9 (*) * (16) 



with 



a I in.-i : x > 



a(x) = { nv n x^ • - ^ » (17) 
: x < 



where the equalities follow methods used in © together with g(A)'s circular- symmetry, and the inequality 
follows from g(\) being monotonically non-increasing together with N'(x) obeying (fl"5l) . 

Define {ai(x)}i e i to be the set of all functions a(x) such that (fl"51) is satisfied. Specifically let's define 
the following two functions: 

• The first function a mg (x) is defined to be piecewise constant over shells defined by consecutive radii 
from the normalized distance series {X^JLq, (i.e. a shell is defined as Sj = {x : Ay_i < x < Xj}). 
The constant value for shell j is selected such that the spectral mass Mj is spread evenly over the 
shell. Formally, this can be expressed as 

[« mg (*)L e5j = v (18) 

• The second function a opt (x) is selected to be etj(x) such that the bound (fTTT) is tightest; thus by 
definition 

j = argmin max aii(x). (19) 

iel x<e 5 Xrt 



4 This restriction holds for many continuous noise channels, such as AWGN. For other channels, it is possible to define similar heuristics. 
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As it turns out a opt (x) is also piecewise constant over shells defined by consecutive radii from the 
normalized distance series, and can be obtained as the solution to a linear program presented in 
appendix [A] 

One subtlety, that can go by unnoticed about a° pt (x), is its dependence on A max . Careful examination 
reveals that a opt (a;) is constructed from those spectral elements whose corresponding distances are 
less than or equal to A max . This is of relevance when optimization of A max is dependent on a opt (a;), 
as is the case in some of the bounds discussed hereafter. This technical subtlety can be overcome by 
construction of a suboptimal version of « opt (x) consisting of more spectral elements than necessary. 
In many cases both the suboptimal and optimal versions coincide. In the remainder of this paper, 
this technicality is ignored and left for the reader's consideration. 



Figure Q] illustrates 9T( 



x), a 



x 



and a° pt (x) for the rectangular lattice 1? . 
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Fig. 1. 0T(x), a' ng (x), and a° vl (x) for the rectangular lattice 



IV. A General ML Decoding Upper Bound 

Many tight ML upper bounds originate from a general bounding technique, developed by Gallager [TT2l . 
Gallager's technique has been utilized extensively in literature 0, [fT3l . 031. Similar forms of the general 
bound, displayed hereafter, have been previously presented in literature (5), |[T5l . Playing a central role 
in our analysis, we present it as a theorem. 

Before proceeding with the theorem, let us define an Additive Circularly-Symmetric Noise (ACSN) 
channel as an additive continuous noise channel, whose noise is isotropically distributed and is a non- 
increasing function of its norm. 

Theorem 3 (General ML Upper Bound). Let A be an n-dimensional lattice, and f\\ z \\ (p) the pdf of an 
ACSN channel 's noise vector 's norm, then the error probability of an ML decoder is upper bounded by 



P e (A)<min( f f\\*\\{p)P2{\p)dp 

\AGA\{0} J ° 



(20) 



+ 



f\\z\\(p)dp 



where the pairwise error probability conditioned on ||z|| = p is defined as 

P 2 (A,p) = Pr(AeBall(z,||z||)|||z||=p) 



(21) 
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where Ball(z, ||z||) z'^ an n-dimensional ball of radius ||z|| centered around z. 

Proof: See appendix El ■ 
We call the first term of (1201) the Union Bound Term (UBT) and the second term the Sphere Bound 

Term (SBT) for obvious reasons. 

In general ((IT)) is difficult to quantify. One method to overcome this, which is limited for analysis 

of specific lattices, is by averaging over the MHS ensemble. New methods for upper bounding (|20l) for 

specific lattices are presented in the next section. 

V. Applications of the General ML Bound 



In this section, the UBT of the ML decoding upper bound (1201) is further bounded using different bound- 
ing methods. The resulting applications vary in purpose, simplicity, and exhibit different performance. We 
present the applications and discuss their differences. 

A. MHS 

Application of the MHS theorem from ©, leads to the random-coding error bound on lattices. Since 
it is based on the MHS ensemble average, the random-coding bound proves the existence of a lattice 
bounded by it, but does not aid in finding such lattice; neither does it provide tools for examining specific 
lattices. 

Theorem 4 (MHS Bound, Theorem 5 of lfT5lO . Let /|| z ||(p) be the pdf of an ACSN channel's noise vector's 
norm, then there exists an n-dimensional lattice A of density (3 for which the error probability of an ML 
decoder is upper bounded by 

/r* poo 
f M (p)p n dp + J ^ f M (p)dp (22) 

with 

r* = (/3K)" 1/n . (23) 

Proof: Set g(X) as 

0(A) = ^ f M (p)P 2 (\p)dp, (24) 
Jo 

noting that it is a bounded function of A and continue to bound the UBT from (l20l) using ©. The 
remainder of the proof is presented in appendix O ■ 

B. DMHS 

Application of the DMHS theorem (fTTI) using a(x) = a opt (x) provides a tool for examining specific 
lattices. The resulting bound is essentially identical to the MHS bound, excluding a scalar multiplier of 
the UBT. It is noteworthy that this is the best a (x) -based bound of this form, since a opt (x) is optimized 
with regards to DMHS. 

Theorem 5 (DMHS Bound). Let a specific n-dimensional lattice A of density (3 be transmitted over an 
ACSN channel with f\\ z \\ (p) the pdf of its noise vector's norm, then the error probability of an ML decoder 
is upper bounded by 

(pr poo \ 

aPV n J f M {p)p n dp + ^ f M {p)dp\ (25) 

with 

a = max a opt (x) (26) 

x<e s -2r 
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where a opt (x) is as defined by (fT9l ). 

Proof: Set g(X) as in (1241) . noting that it is bounded by A max = 2r. The remainder is identical to the 
proof of theorem |H replacing f3 with a (3. ■ 
Optimization of r can be performed in the following manner. Since a opt (x) is a monotonically non- 
increasing function of x, optimization of r is possible using an iterative numerical algorithm. In the 
first iteration, set r = (/3V n )~ 1 / n and calculate a according to (l26t . In each additional iteration, set 
r = (ctfiVn)' 1 ! 1 ' 1 and recalculate a. The algorithm is terminated on the first iteration when a is unchanged. 



C. eDMHS 

Rather than maximizing the UBT using a single scalar factor (as was done in the DMHS), the eDMHS 
splits up the UBT integral to several regions with boundaries defined by {A-,}°^ . Maximization of each 
resulting region by its own scalar, results in much tighter, yet more complex bound. This is typically 
preferred for error bounding in finite dimension lattices. For asymptotical analysis, the eDMHS would 
typically require an increasing (with the dimension) number of spectral elements, while the DMHS would 
still only require one, making it the favorable. We choose a(x) = a opt (x), rather than optimizing a(x) for 
this case. The reason is that although this bound is tighter for the finite dimension case, it is considerably 
more complex than a competing bound (for the finite dimension case) presented in the next subsection. 
The main motivation for showing this bound is as a generalization of DMHS. 

Theorem 6 (eDMHS Bound). Let a specific n- dimensional lattice Ao of density (3 be transmitted over 
an ACSN channel with f\\ z \\(p) the pdf of its noise vector's norm, then the error probability of an ML 
decoder is upper bounded by 

f\\z\\(p) h j(p)dp 

V2 




fhwWP) ( 27 ) 



with 

otj = maxa opt (s) = a opt (e%) (28) 



hj(p) 



[ a{xGBall(z, ||z||)| ||z|| = p}rfx (29) 
J Si 

Sj = {x :x> 0, e%_i < x < e%} (30) 

Sj = {x : x G E n , Aj-_i < ||x|| < Xj} (31) 

where a opt (x) is as defined in (fT9l) , cr{x G Ball(z, ||z||)} is the characteristic function of Ball(z, ||z||), 
{Aj}°^ is the previously defined distance series of A , and M is the maximal index j such that Xj < 2r. 

Proof: We set g(X) as in (|24|) noting that it is bounded by A max = 2r. We continue by remembering 
that a opt (x) is piecewise constant in the shells Sj, and therefore (fT2)) collapses to 

M „ 

gW<PY, a i / 9^ d *- (32) 

AeA /{0} 3=0 Js i 

For the remainder we continue in a similar manner to the proof of theorem |4] by upper bounding the UBT 
from (|20l ) using (1321 . See appendix ID1 ■ 
A geometrical interpretation of hj(p) is presented in figure[2l 
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D. Sphere Upper Bound (SUB) 

Another bound involving several elements of the spectrum can be constructed by directly bounding 
the conditional pairwise error probability P 2 (x,p). This bound is considerably more complex than the 
DMHS, but is potentially much tighter for the finite dimension case. When the channel is restricted to 
AWGN, the resulting bound is similar in performance to the SUB of H, hence the name. The bounding 
technique for P 2 (x, p), presented in the following lemma, is based on lfT6ll . 



Lemma 2 (Appendix D of aim ). Let ~x.be a vector point in n-space with norm ||x|| < 2p, z an isotropically 
distributed n- dimensional random vector, and p a real number then 




Pr(xeBall(z.||z||)| ||z|| = ,>) < ( 1 - ( ^ ) j (33) 

Proof: See appendix [0 ■ 
The above lemma leads directly to the following theorem. 

Theorem 7 (SUB). Let a specific n-dimensional lattice Ao of density (3 be transmitted over an ACSN 
channel with /|| z ||(p) the pdf of its noise vector's norm, then the error probability of an ML decoder is 
upper bounded by 

P e (A ) < min 





(34) 

where and {Mj}'jL 1 are the previously defined distance series and spectrum of A respectively, 

and M is the maximal index j such that Xj < 2r. 

Proof: Bound the UBT of (HI directly, using <|33]). See appendix E ■ 
Optimization of r can be performed in the following manner. We begin by analyzing the function 

/ 2\ "2 1 

f(p) = ( 1 — (fy) J ■ It is simple to verify that this function is positive and strictly increasing in the 

domain {p : p > Xj/2}. Since the UBT of (l34l) is a positive sum of such functions, it is positive and 
monotonically nondecreasing in r. Since additionally the SBT of (l34l) is always positive, it suffices to 
search for the optimal r between Ax/2 and an r max , where r max is defined such that the UBT is greater 
than or equal to 1. 
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We continue by calculating M max that corresponds to r max . By definition, each selection of M corre- 
sponds to a domain {r : Am < 2r < Am+i}- Instead of searching for the optimal r over the whole domain 
{r : Ai/2 < r < r max }, we search over all sub-domains corresponding to 1 < M < M max . When M is 
constant, independent of r, the optimal r is found by equating the differential of (|34|) to 0, in the domain 
{r : Am < 2r < Am+i}- Equating the differential to results in the following condition on r: 

J -1 = 0- <35) 

The function on the left of the above condition is monotonically nondecreasing (as shown previously), so 
an optimal r exists in {r : Am < 2r < Am+i} iff its values on the domain edges are of opposite signs. If 
no optimum exists in the domain, it could exist on one of the domain edges. 

The optimization algorithm proceeds as follows: set M = 1 and examine the differential at the domain 
edges Am/2 and Am+i/2. If the edges are of opposite signs, find the exact r that zeros the differential 
in the domain and store it as r^f. Otherwise set r^f = 0, advance M by 1 and repeat the process. The 
process is terminated at M = M max . When done, evaluate (1341 ) at r = r\ cl . . . , r^f , Ai/2, . . . , Am 1mx /2 and 
select the r that minimizes. 

We conclude this section with an example presented in figure [3l which illustrates the effectiveness of 
the new bounds for finite dimension. The error probability of the Leecrjf] lattice A 2 4 is upper bounded 
by DMHS <g5) and SUB d34]). The ordinary Union Bound (UB), the MHS bound for dimension 24 (1221) . 
and the Sphere Lower Bound (SLB) of [fTTl are added for reference. The spectral data is taken from [fT8l . 
The bounds are calculated for an AWGN channel with noise variance a 2 . 




VNR=e~ 25 /o 2 



Fig. 3. A comparison of DMHS and SUB for the Leech lattice. The UB, MHS and SLB are added for reference. The graph shows the 
error probability as a function of the Volume-to-Noise Ratio (VNR) for rates S* < 5 < 5 cr , 



VI. Error Exponents for the AWGN Channel 

The previous section presented three new spectrum based upper bounds for specific lattices. As stated 
there, only the DMHS bound seems suitable for asymptotical analysis. This section uses the DMHS 
bound to construct an error-exponent for specific lattices over the AWGN channel. Although this case 
is of prominent importance, it is important to keep in mind that it is only one extension of the DMHS 
bound. In general, DMHS bound extensions are applicable wherever MHS bound extensions are. 

5 The Leech lattice is the densest lattice packing in 24 dimensions. 
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When the channel is AWGN with noise variance a 2 , the upper error bound on the ML decoding of 
a "good'H lattice from the MHS ensemble (1221) can be expressed in the following exponential form (3), 

ma 

P e (A) < e -n(MS)+o(i)) ( 36) 

with 

(<f*-<f)+logf, 5<5 cr 
£ r (5) = <J gt^ z |g z fl z l ) 5cr < 5<( 5* (37) 
0, " 5 > 5* 

5* = \ log - (38) 
5 cr = ilog-^ (39) 

where o(l) goes to zero asymptotically with n. 

This error-exponent can be directly deduced from the MHS bound (l22l) . By applying similar methods 
to the DMHS bound (1251) . it is possible to construct an error-exponent for a specific lattice sequence based 
on its distance- spectrum. 

Theorem 8 (Non-Random Coding Error Exponent). Let A [n] be a specific lattice sequence transmitted 
over an AWGN channel with noise variance a 2 , then the error probability of an ML decoder is upper 
bounded by 

Pe(AoN) < e -n(M6+"ln])+o(l)) (40) 



A 1 



with 

v[n] = — logafn]. (41) 

where [n] indicates the n'th element of the sequence. 

Proof: It follows from the proof of theorem \5\ that replacing f3 with /3a [n] there, is equivalent to 
replacing 5 with S + u[n] here. ■ 

Clearly (l40l) can be used to determine the exponential decay of the error probability of a specific lattice 
sequence. This leads to the following corollary. 

Corollary 2 (Gap to Capacity 1). A lattice sequence for which 

lim - logafra] = (42) 

rwoo n 

achieves the unrestricted channel error-exponent. 



Proof: Follows immediately from (1401) and the definition of u[n\. ■ 
When certain stronger conditions on the lattice apply, the following simpler corollary can be used. 

Corollary 3 (Gap to Capacity 2). A lattice sequence for which a mg (x) (per dimension n) is monotonically 
non-increasing in x and 

lim 1 log )=0 (43) 

n->mfi \e n V n {Ai[n\) n J 

achieves the unrestricted channel error-exponent. 

Proof: When a rng (x) is monotonically non-increasing then a opt (x) = ct rng (a;) and a[n] = ^tt^t^ • 

■ 

6 Define a "good" lattice, from an ensemble, as one that is upper bounded by the ensemble's average. 
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Let us examine corollary [3j Assume a sequence of lattices with minimum distance Xi[n] = e~ 5 Vn 
is available (see [fT9ll for proof of existence). Plugging this back into (|43l ) leads to the following necessary 
condition for the bound (|40l) to achieve the unrestricted channel error-exponent 

lim - log (MW) = 0, (44) 

n-J-oo Tl 

whether the monotonicity condition applies or not. Although this only amounts to conditioning on the 
bound (and not on the lattice sequence), condition (1441) gives an insight to the close relationship between 
the spectral distances and their enumeration. We conjecture that at rates close to capacity, the bound is 
tight, leading to (|43l) being a necessary condition on the lattice sequence itself. 

We conclude this section with an illustration of the exponential decay series u[n) of the first three 
lattices of the Barnes-Wall lattice sequence BW4 = -D4, BW& = E 8 , and BW\§ = Ai 6 . Unfortunately 
the distance-spectrum for BW n is generally unknown, preventing asymptotical analysis. Nonetheless an 
interpolation for dimension 4 to 16 is presented in figure 51 The spectral data is taken from |fT8l . 

0.25 
0.2 
0.15 

IE 

> 

0.1 
0.05 



"4 6 8 10 12 14 16 

n 

Fig. 4. The exponential decay series u[n] for the lattice sequence BW„, calculated for dimensions 4, 8, and 16 and interpolated in-between. 

Examination of figure |4] shows that the upper bound on the gap to capacity decreases with the increase 
in n, at least for the first three lattices in the sequence. Although full spectral data is not available for the 
remainder of the sequence, the minimal distance and its enumeration are known analytically. Assuming 
momentarily that the condition on a mg (x) as presented in corollary [3] holds, we can try to examine 

n ( e^v^l^ln])" ) as an u PP er bound on the gap to capacity. This examination is illustrated in figure 
|5l where the first three lattices coincide with the previous results in figure 0] Clearly, these results are 
a lower bound on an upper bound and cannot indicate one way or the other. Nonetheless, it seems that 
the results are consistent with the well-known coding performance of Barnes-Wall lattices with increasing 
dimension. 

Appendix 

A. Linear Program for Finding a opt (x) 

We begin by restating the conditions describing a opt (x). 

a °*(x) = aj {x) 

j = argmin max «j(x) (45) 
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log 2 n 

Fig. 5. i log ( ^rffiy^nfyi J f° r the first few dimensions of the lattice sequence BW n . Coincides with i log a[n] for the first three 
lattices. 



where a(x) is expressed by 



N'(x) 



:x^0 



a{x) A ^ - ^ u (46) 



and N'(x) satisfies the inequality 

pr 

N(x)dx< I N'(x)dx Vr G (0,e 5 A max ]. (47) 



Recall that N'(x) is a smoothed version of the continuous spectrum N(x) such that spectral mass is allowed 
to shift down from higher to lower distances. As is apparent from the conditions above, finding a opt (x) 
is equivalent to finding N'(x) such that (|45T) is optimized. It is immediately clear from (|45l) that since we 
are searching for the function that minimizes its maximum value over the domain {x : < x < e 5 \ max }, 
we can only benefit from limiting N'(x) to consist of only the spectral elements {A/}}^ where M is the 
maximal index such that Xj < A max - This limiting necessarily does not increase the max value over the 
domain {x : < x < e^Amax} and is consistent with (|47T) for r = A max . Clearly, optimization of (1451) is 
achieved by an a(x) that is as constant as possible over the domain {x : < x < e^A^x}, as is typical for 
a minmax optimization. We say "as constant as possible" since (l47l) imposes constraints on our ability 
to find N'(x) such that a(x) is constant over the whole domain {x : < x < e^Amax}. Given (|47l) . the 
optimum a(x) is piecewise constant over the sub-domains {x : Aj-_i < x < Xj} with 1 < j < M, and is 
zero in the sub-domain {x : \m < x < A max }. With the piecewise nature of a(x) established, and due to 
d46l) 

[ 3 N'{x) = I 1 a{x)nV n x n - l dx 
= ctj nV n x n ~ 1 dx 

= a 3 V n (\]-X r ;_^ (48) 

where aj is the constant value of a(x) over sub-domain {x : Aj__i < x < It is quite clear from (l48~l) . 
that in our case N'(x) can be completely described by tj = j~ 3 N'(x) for 1 < j < M. 
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We now present the general formulation necessary to construct the LP for optimizing a(x). Take each 
spectral element Mj and subdivide it into j parts {A/},»}| =1 such that A/},, > and Mj = YH=\ Nj,i- Mj t i can 
be viewed as the contribution of the j'th spectral element Mj to the i'th sub-domain {x : Aj_x < x < Aj} 
where 1 < % < j. Defining to be the total spectral mass associated with sub-domain {x : Aj_i < x < Aj}, 
we can summarize the description of the LP variables A/},j as follows: 

> o 



E-^ 

i=l 

M 

E^ 



Mi 



V i<j<M 

V 1 < j < M 



ei V 1 < j < M 

An illustration of the LP variables for M = 4 is presented in the following table: 



(49) 
(50) 

(51) 





M 


M 2 


M 3 


M 


ei 




M 2jl 


M 3 ,i 








Af 2 ,2 






^3 






A^3,3 




C4 








A/4,4 



The LP variables set {{A/}^}^}*!! completely defines N'(x), such that its optimum completely defines 
° pt (x). Finally the LP is defined as follows: 



minimize a 

M jti > 



E^=^ 



V 
V 



i<j<M 
1 < J < M 



i=l 



E 



M M„ 



3=1 



3.1 



V n A!- 



< a 



V 1 < j < M. 



(52) 
(53) 

(54) 



It is easy to show that the LP optimization landscape is bounded, by bounding each variable individually 

< Mj,i < Mj V i<j<M. (55) 
It is simple to show that there exists a solution to the LP constraints. One trivial solution is 



M 3 




: 1 = j 



V i<j<M. 



(56) 



The LP can be solved by a water-filling algorithm which is linear in the number of considered spectral 
elements M . At the j'th step the water (a.k.a. spectral mass Mj) is poured into the j'th bucket (a.k.a. sub- 
domain {x : Aj_i < x < \j}), until it overflows (a.k.a. exceeds the previous maximum), at which point 
it is permitted to pour over lower sub-domains. A graphical representation of the algorithm is presented 
in figure |6l The optimality of the water-filling solution can be verified by induction. 

B. Proof of Theorem [£] 

We use Herzberg and Poltyrev's upper sphere bound definition for AWGN [6] to introduce the general 
technique in the context of lattices over ACSN channels. Consider an n-dimensional lattice A transmitted 
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Fig. 6. A graphical representation of the water-filling solution. 



over an ACSN channel. In each transmission a lattice point A is transmitted and y = A + z is received. 
The additive noise z corrupts the transmitted lattice point. An ML decoder decodes correctly so long as 
z is contained in the voronoi cell surrounding A. The probability of decoding error can be expressed as 

P e (A) = Pr(e| || z || < r) Pr(||z|| < r) 

+ Pr(e| ||z|| > r) Pr(||z|| > r) (57) 



where e denotes a decoding error event, and r is a real positive parameter. Since Pr(-) < 1, equation (|57l) 
can be trivially upper bounded as 

P e (A) < Pr(e| ]| zs || < r) + Pr(||z|| > r). (58) 



Assuming without loss of generality that A = was transmitted, the first term of (1581) can be upper 
bounded as 

Pr(e| ||z|| < r) = Pr J (ex| ||z|| < r) 

\AeA/{0} 

< Pr(e A | ||z|| < r) (59) 

AeA/{0} 

where e\ is a pairwize error resulting from decoding A ^ 0, and the inequality is due to the union bound. 
The bound after substitution of (|59l) in (1581) is given by 

P e (A) < Yl Pr (e A | ||z|| < r) + Pr(||z|| > r). (60) 

AeA/{0} 

Consequently, the tightest bound is found by optimization over r, 

P e (A)<min ^ Pr(e A | ||z|| < r) + Pr(||z|| > r) | . (61) 

\AGA/{0} 
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Following techniques from Ingber et al. lfT5ll , the inner probability term of the first term on the right- 
hand-side of (|60l ) can be expressed as^| 

Pr(eA| ||z|| < r) 

= / /|]z[iO)Pr(e x | ||z|| = p)dp 
Jo 

= [ /| W |(p)Pr(||z-A||<||z|||||z||=p)dp 
Jo 

= [ / N |(p)Pr(AeBaU(z,||z||)|||z||=p)dp (62) 
Jo 

where f\\ z \\(p) is the pdf of the noise vector's norm, and Ball(z, ||z||) is an n-dimensional ball of radius 
||z|| centered around z. Using the same notations, the second term on the right-hand-side of (l60l) can be 
expressed as 

poo 

Pr(||z|| >r) = J f M {p)dp. (63) 
Plugging (l62l) and (l63l) into (|6TT) completes the proof. 

C. Proof of Theorem |4] 

This proof is brought for completeness and reading convenience, although it has been seen, in similar 
form, elsewhere, previously. 
Using inequality ©, 



f /||z||(p)Pr(A G Ball(z, ||z||)| ||z|| = p)dp 
AeA/{0} Jo 

<P f I /ii-ii(p) Pr ( x e Ball ( z , IN|)| ||z|| = P )dpd* 
Ju n Jo 

= P I /|M|(P) / Pr (x e Ball(z, ||z||)| ||z|| = p)d*dp 
Jo Ju n 

= pf /n(p) / E(cr{x G Ball(z, ||z||)}| ||z|| = p)dxdp 
Jo Ju n 



p dp 



P f /|H|(p)E( / a{xGBall(z,||z||)}rfx ||z| 

Z 3 / /|M|(p) / a { x G Ball ( z , ll z ll)l ll z ll = p}d*dp 
Jo Jn n 

0V n f f M {p)p n dp (64) 



o 



where a{x G Ball(z, ||z||)} is the characteristic function of Ball(z, ||z||). The optimal value r* of r is 
found by differentiating (1221) . 

PV n f M (r)r n = f M (r). (65) 

7 The integral in l !62l l can be further constrained to [||A|| /2, r] since Pr(A G Ball(z, ||z||)j z = p) — for p < ||A|| /2. 
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D. Proof of Theorem [6] 
Using inequality (|32|) . 



o 



£ / /||z||(p)Pr(AGBall(z,||z||)|||z||=p)rfp 

AGA /{0}' " 
M 

< 

3=0 



A3 « nj" 

P^atj / /|| z ||(p) Pr(x G Ball(z, ||z||)| ||z|| = p)dpdx. 

„-_n ^5, 7o 



(66) 



Continue with the term inside the sum (|66l) . 



/ / /,| z ||(p)Pr(xGBall(z,||z||)|||z||=p)rfpdx 

/]| z ]|(p) / Pr(x G Ball(z, ||z||)| ||z|| = p)dxdp 
f M (p) [ E(a{xGBaU(z,||z||)}|||z||=p)dxdp 

/ h (p)eQ a{x G Ball(z, ||z||)}dx 

/ ^{x G Ball(z, ||z||)| ||z|| = p}dxdp 

fM(p) h j(p) d P 
f\\z\\(p)hj(p)dp 



|z|| = p ] dp 



(67) 



V2 



where a{x G Ball(z, ||z||)} is the characteristic function of Ball(z, ||z||), and the last two equalities follow 
from the definition of hj(p). Plugging (|67l) back into (|66l) completes the proof. 



E. Proof of Lemma [2] 

Clearly when ||x|| > 2p, the probability Pr(x G Ball(z, ||z||)| ||z|| = p) = 0. 
Begin by defining the i.i.d. Gaussian vector 

u~A/" n (0,l). 

The vector z conditioned such that ||z|| = p can be expressed as 



-u 



(68) 



(69) 



with no loss of generality. Since u is isotropically distributed, the vector x can be set to 

x= [||x||, (),••• ,0] (70) 
again with no loss of generality. Continue by noting that x G Ball(z, ||z||) is equivalent to ||z — x|| < p, 
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thus 



Pr(x e Ball(z, ||z||)| ||z|| = p) 
= Pr(||z-x|| <p) 



= Pr 
= Pr 



P 



||z| 



[zi, 



|x||,0,... ,0] 



< 



P 



Pr U - 



|x|| ||z| 
P 

Ixll llzl 



j z 2, 
2 



< Z 



Pr 2i > 



P 

Ixll ||z| 



2p 



= -PrU> 



|x|| /2p) s 



l-N|/2p) : 



E 

i=2 



Substitute £ = (||x|| /2p) 2 for easier readability, thus 



1 

2 



^Pr ( %\ > — - \\z 2 :n\\ 2 



E^-Pr( z\> 



1-t 



Z 2: . 



Z2:>, 



t 



< E exp 



1-t 
1 t 



\Z2:n\ 



21 -t 



Z 2 :: 



= J i exp ^-1^— Ik2:nf^ fAf^(0,l)(z 2 :n)dz 2 :n 



= (!"*) 
= (l-t) 



n-l 
2 



n-l 
2 



/aT™- 1 ^,!-*) (Z2:n)dz2:n 



(71) 



(72) 



where the inequality results from a fairly crude upper bound on the Gaussian error function = 
Pr(A/"(0, 1) > x) < exp(— x 2 /2). Undoing the substitution of i completes the proof. 
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F. Proof of Theorem [7] 



f f\w(p) Pr(A E Ball(z, ||z||)| ||z|| = p)dp 

AeA /{0} Jo 

f f\HiP) Pr ( A e Ball (z, Nl)| |M| = P)dp 



AeA /{0} 1 



^ E / /H(p)(i-(y) ) ^ 



AeA /{o} >/ ll A ll/ 2 



n-1 
2\ ~~2~~ 



n-1 

M „ r / / x \ 2\ ~2 



= E^i>HW^-(8) * C73) 

where the first equality is due to Pr(A e Ball(z, ||z||)| ||z|| = p) = for p < ||A|| /2, the inequality is due 
to (1331) . and the last equality is formed by aggregation of identical distance summation terms. 
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